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Abstract 


A  method  of  collecting  ensembles  for  conditional  averaging  is  presented  that 
uses  data  collected  from  a  plane  mixing  layer.  The  correlation  coefficient  of 
a  sine  function  and  a  reference  signal  were  used  to  determine  the  adequacy 
of  blocks  of  the  reference  signal  and  the  phase  alignment  for  the  data. 
Selection  of  the  sine  function  period  and  a  correlation  coefficient  threshold 
are  discussed.  Also  examined  are  the  effects  of  the  period  and  threshold  level 
on  the  number  of  ensembles  captured  for  inclusion  for  conditional  averag¬ 
ing.  Both  the  selection  of  threshold  correlation  coefficient  and  the  choice  of 
data  set  length  determine  bandpass  filtering;  this  study  examined  the  effect 
of  these  choices  on  the  ability  to  discriminate  reference  velocity  ensembles 
for  conditional  averaging. 
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1.  Introduction 

One  powerful  tool  for  investigating  coherent  structures  in  a  turbulent  fluid 
flow  is  a  conditional  average  resulting  from  the  collection  of  several  condi¬ 
tional  samples  that  are  subjected  to  some  averaging  technique.  The  pur¬ 
pose  is  to  determine,  given  enough  events,  what  properties  of  turbulent 
fluid  flow  are  consistent,  and  to  possibly  identify  coherent  structures  in  an 
otherwise  chaotic  flow.  To  perform  a  conditional  average,  one  must  iden¬ 
tify  an  event  in  the  flow,  and  some  feature  about  the  event  with  which  to 
align  the  data  to  be  collected.  An  example  of  an  event  is  the  passage  of  a 
shedding  vortex  off  a  bluff  body.  In  this  example,  the  data  surrounding  the 
event  would  be  aligned  with  the  point  in  time  at  which  the  velocity  is  at  a 
maximum  when  a  vortex  passes.  Once  an  event  is  identified  and  a  means 
of  alignment  provided,  several  separate  events  can  be  measured  and 
pooled  for  conditional  averaging  and  evaluation.  Subramanian  et  al1  de¬ 
scribe  several  techniques  for  detecting  and  averaging  data  based  on 
discriminable  events.  Davies2  used  peak  velocity  as  a  means  to  detect  the 
passage  of  shedding  vortices  in  a  Von  Karman  street,  which  have  a 
strongly  periodic  flow.  In  the  case  of  the  turbulent  mixing  layer,  another 
strongly  periodic  flow,  it  is  obvious  from  inspection  that  the  passage  of  a 
roller  is  an  excellent  candidate  for  an  event.  In  fact,  Hussain  and  Zaman3 
showed  that  a  single  sensor  probe,  placed  several  momentum  thicknesses 
from  the  mixing  layer  centerline,  resulted  in  a  streamwise  velocity  time  se¬ 
ries  that  exhibited  clear,  periodic  events.  As  in  the  mixing  layer  depicted  in 
figure  1,  the  single  hot-wire  probe  that  is  located  at  the  high-speed  edge  of 
the  roller  senses  a  decrease  in  local  velocity  induced  by  the  retarding  mo¬ 
tion  of  the  passing  roller.  When  the  roller  passes,  the  local  streamwise  ve¬ 
locity  increases  above  that  of  the  free  stream.  The  velocity  oscillates,  and 
appears  nearly  sinusoidal  with  the  passage  of  time. 


Figure  1.  Passage  of 
roller  with  respect  to 
reference  probe. 
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1C.  S.  Subramanian,  S.  Rajagopalan,  R.  A.  Antonia,  and  A.  ].  Chambers,  "Comparison  of  conditional  sampling  and 
averaging  techniques  in  a  turbulent  boundary  layer,"  J.  Fluid  Mech.  123  (1982),  335-362. 

2M.  E.  Davies,  "A  comparison  of  the  wake  structure  of  a  stationary  and  oscillating  bluff  body,  using  a  conditional 
averaging  technique,"  J.  Fluid  Mech.  75  (1976),  209-231. 

3A.K.F.M.  Hussain  and  K.B.M.Q.  Zaman,  "An  experimental  study  of  organized  motion  in  the  turbulent  plane  mix¬ 
ing  layer,"  J.  Fluid  Mech.  159  (1985),  85-104. 
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Figure  2.  Fast  Fourier 
transform  (FFT)  of 
reference  velocity 
power  spectrum, 
experiment  1. 


We  can  take  advantage  of  this  sinusoid-like  signal  to  detect  the  passage  of 
the  roller,  and  collect  the  data  of  interest  surrounding  the  event.  We  would 
align  the  data  merely  by  identifying  some  points  in  the  signal,  such  as  the 
maxima  or  minima.  Unfortunately,  the  velocity  data  are  not  exactly  a  sine 
function.  In  fact,  as  seen  in  the  fast  Fourier  transform  (FFT)  of  figure  2, 
there  is  a  range  of  frequencies  associated  with  roller  passage.  An  alterna¬ 
tive  method  for  filtering  the  reference  velocity  and  aligning  the  data  is 
needed. 


Roller  passage  frequency  band 
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2.  Event  Detection 

One  method  to  detect  the  passage  of  a  particular  type  of  roller  is  to  use  a 
correlation  of  the  reference  velocity  with  a  sine  function,  as  depicted  in 
figure  3a.  The  correlation  of  the  data  with  a  test  sine  function  results  in  a 
correlation  coefficient: 


\ref~ 


UtUref 


V“F* 


ref 


(1) 


where  rt  rejis  the  correlation  coefficient,  ut  is  the  test  sine  function  sin(2rojf), 
and  uref  is  the  reference  velocity.  The  test  function  is  correlated  with  a  seg¬ 
ment  of  data  of  equal  length,  T,  and  the  correlation  coefficient  is  calculated. 
The  test  signal  is  then  correlated  with  the  data  that  have  been  shifted  by 
some  At,  which  is  very  small.  If  the  data  are  made  up  of  a  series  of  data 
points,  such  as  digital  data,  then  the  shift  consists  of  one  point.  Thus,  the 
entire  original  set  (except  the  very  first  point),  plus  a  new  point  at  the  end 
of  the  set,  is  correlated  with  the  test  function,  and  a  new  correlation  coeffi¬ 
cient  is  calculated.  Eventually,  as  the  reference  velocity  signal  is  advanced 
and  correlated,  a  time  history  of  rt  ref  is  created,  as  in  figure  3b.  This  his¬ 
tory  of  rt  rej  enables  us  to  accomplish  two  objectives:  the  detection  of  the 


Figure  3.  Reference 
roller  velocities: 

(a)  reference  velocity 
fluctuations  compared 
to  sine  function 
segment  and 

(b)  correlation 
coefficient  with  signal 
segment. 
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desired  segment  of  reference  velocity,  and  the  alignment  of  the  data 
sample,  as  the  correlation  coefficient  will  reach  a  peak  and  then  reduce  as 
the  roller  passes  on.  The  data  would  then  be  collected  for  inclusion  into  the 
conditional  averaging. 

The  likelihood  of  encountering  a  correlation  of  r  =  1.0  is  very  small.  In  the 
present  study,  and  as  evident  from  the  data  sample  in  figure  3b,  this  never 
occurred.  A  threshold  for  r  had  to  be  established  to  provide  enough  condi¬ 
tional  samples  of  the  reference  velocity  to  create  a  good  statistical  popula¬ 
tion.  Since  there  was  a  band  of  frequencies  in  the  reference  velocities,  the 
optimal  frequency  for  the  sine  function  also  had  to  be  picked. 


3.  Determination  of  Correct  Test  Frequency 

The  object  in  selecting  a  test  frequency  is  to  capture  as  many  realistic  en¬ 
sembles  as  possible  from  a  data  set.  The  test  frequency  may  initially  be  de¬ 
termined  from  the  power  spectrum  plot  of  the  reference  velocity  from  fig¬ 
ure  4.  The  frequencies  ranged  from  2  to  10  Hz,  with  a  peak  at  about  6  Hz. 
This  indicates  that  the  energy  distribution  has  a  concentration  in  a  periodic 
of  that  frequency.  To  determine  the  optimal  frequency  for  the  test  signal,  I 
performed  a  series  of  tests.  The  reference  velocity  data  were  correlated 
with  a  variety  of  sine  functions  with  different  values  of  a.  The  threshold 
correlation  coefficient  r ^  was  held  to  0.5.  This  threshold  coefficient  seemed 
appropriate.  The  sample  spanned  4000  points,  or  0.4  s  of  data.  This  sample 
size  was  used  since  it  was  desired  to  capture  at  least  one  complete  roller, 
plus  most  of  the  data  associated  with  the  rollers  preceding  and  following 
the  central  roller.  As  each  test  frequency  was  used,  the  total  number  of  en¬ 
sembles  captured,  N,  was  recorded.  Figure  5  shows  the  plot  of  N  to  a  for 
the  first  experiment  reference  velocity.  The  peak  number  of  captures  ap¬ 
peared  to  be  at  about  a  =  8.5,  which  was  different  from  the  peak  frequency 
in  figure  4. 


Test  frequency,  a  (Hz) 


Figure  6.  Effect  of 
normalization  on  test 
frequency. 


The  large  number  of  ensembles  accumulated  at  a  frequency  other  than  at 
the  power  spectrum  peak  is  due  to  the  value  used  for  the  correlation 
threshold  rth.  The  value  of  the  threshold  would  have  a  direct  impact  on  the 
reference  velocity  bandwidth.  A  rudimentary  method  of  eliminating  the 
impact  of  rth  would  be  to  normalize  N  by  the  total  number  of  expected 
roller  passages  in  the  data  set,  if  uref  behaved  as  ut.  The  expected  number  of 
ensembles  for  urer  would  be  ciTset,  where  Tset  is  the  time  length  of  the  entire 
data  set.  Normalization  would  be  as 


Figure  6  shows  the  effect  of  this  normalization.  The  peak  is  shifted  back  to 
the  value  indicated  by  the  power  spectrum.  This  indicates  that  looking  at 
N  alone  will  not  enable  selection  of  the  appropriate  test  frequency.  Care 
must  be  given  to  selecting  the  proper  threshold  correlation  coefficient  rth. 
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4.  Selection  of  Threshold  Correlation  Coefficient 

Selection  of  the  threshold  correlation  coefficient  rth  affects  the  number  of 
roller  passages  that  will  be  included  in  the  conditional  average.  Raising  or 
lowering  this  value  is,  in  effect,  the  same  as  changing  the  bandwidth  of  a 
band-pass  filter.  The  bandwidth  would  be  centered  about  the  test 
frequency  a,  and  rollers  with  associated  reference  velocity  that  have 
passed  the  frequency  bandwidth  would  be  included  for  conditional  aver¬ 
aging.  To  estimate  the  bandwidth  as  a  function  of  rth,  I  compared  the  test 
signal  with  another  function,  sin(2 nbt).  The  variance  of  a  sine  function  is 
known:  (0.5)“1/2.  Therefore,  the  correlation  function  for  the  continuous 
functions  is  modified  from  equation  (1)  to 

Y  f  sin(27Mf)  sin(27ibf)  dt 

r(a,b,T)  =  — - Q-jj -  •  (3) 


This  evaluates  to 


,  .  -  2/l  sin [2nT(a -b)}  1  sm[2*T(a  +  b)]\  (  ) 

r(a,b,T)  =j\2 - ^ - +  2  ^Tb  /  '  W 

If  we  normalize  b  and  the  sample  length  T  by  the  test  signal  frequency, 
which  is  a  constant  a,  we  obtain  the  response  function  in  normalized  terms: 


Q  =  aT  , 


(R cs\  -  1  [sin  [In&'.  l  -  Pj)  sin  [2tcO(1  + /?)]) 

r{P,)~Q\  l-p  1  +  P  /' 

Equation  (5)  can  be  used  to  determine  the  bandwidth  of  accepted  normal¬ 
ized  frequencies  about  p  =  1  for  a  given  r^.  It  also  shows  that  bandwidth  is 
not  solely  dependent  on  frequency,  but  also  on  the  size  of  the  sample  space 
T.  Figure  7  shows  the  contour  plot  of  r  using  the  normalized  variables.  The 
reference  frequency  P  is  varied  along  the  vertical,  while  0  is  varied  along 
the  horizontal.  The  contours  are  values  of  r. 

Use  of  equation  (5)  helps  to  select  the  sample  length  and  the  associated 
threshold  correlation  coefficient.  Figure  8  shows  an  example  of  bandwidth 
from  given  values  of  0  and  r. 

With  the  threshold  correlation  coefficient  set  to  0.6,  the  range  of  frequen¬ 
cies  passed  is  P  =  0.899  to  1.096.  Figure  9  shows  the  A p  bandpass  range 
with  respect  to  the  threshold  correlation  coefficient. 
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We  still  need  to  determine  the  appropriate  test  frequency,  threshold  corre¬ 
lation  coefficient,  and  sample  space.  In  the  example  of  figures  8  and  9,  a 
sample  space  of  4000  points  was  used  (0.4  s  of  data),  with  the  test  fre¬ 
quency  at  6.7  Hz,  which  is  the  peak  value  derived  from  figure  6.  The  band¬ 
width  at  rth  =  0.6  is  A/=  1.32  Hz,  or  a  frequency  range  of  6.02  to  7.34  Hz.  In 
figure  4,  the  center  of  the  region  where  the  roller  passages  contain  the 
greatest  spectral  energy  is  at  about  6  Hz.  This  high-energy  region  ranges 
from  approximately  4  to  7  Hz.  Assuming  the  time  scale  is  smaller,  fixed  at 
T  =  0.2  s,  we  can  determine  the  approximate  Af  from  the  spectral  plot,  and 
evaluate  the  lowest  value  of  rtj,.  Assuming  that  the  desired  range  would  be 
3  Hz,  centered  on  6  Hz,  A/3  is  then  0.5,  or  /3  =  1.0  ±  0.25,  and  0  =  1.2.  From 
figure  7(a)  or  equation  (5),  it  is  found  that  r ^  should  not  be  lower  than 
0.497.  This  threshold  value  would  result  in  a  bandpass  of  A/3  =  0.25,  so  that 
the  correlation-based  conditional  averaging  would  conservatively  capture 
the  desired  information,  while  filtering  out  information  with  spectral  en¬ 
ergy  outside  the  desired  bounds. 
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5.  Conclusion 


The  method  of  roller  passage  detection  used  by  Hussain  and  Zaman3  is 
frequent  in  the  literature.  Few  of  the  methods  described  include  any  pro¬ 
cess  to  discriminate  the  types  of  rollers  that  pass  (Subramanian  et  al1 *).  In 
the  self-similar  region  of  die  mixing  layer,  rollers  sometimes  merge  or  pair. 
The  evolution  of  roller  pairing  involves  a  period  in  which  the  dissipation 
rate  and  turbulent  kinetic  energy  levels  are  higher  than  during  the  period 
in  which  a  roller  is  simply  a  pure  roller.  Inclusion  of  these  types  of  activi¬ 
ties  into  the  conditional  average  would  skew  the  distributions  away  from 
that  of  a  pure  roller.  Use  of  the  correlation  coefficient  effectively  filters  out 
roller  pairings  and  rollers  that  may  not  conform  to  what  may  be  consid¬ 
ered  normal.  Conversely,  the  correlation  coefficient  can  be  used  to  exclude 
pure  rollers,  so  that  one  can  concentrate  on  those  that  are  about  to  pair. 
More  work  needs  to  be  done  with  this  concept. 


lC.  S.  Subramanian,  S.  Rajagopalan,  R.  A.  Antonia,  and  A.  J.  Chambers,  "Comparison  of  conditional  sampling  and 

averaging  techniques  in  a  turbulent  boundary  layer,"  J.  Fluid  Mech.  123  (1982),  335-362. 

3A.K.F.M.  Hussain  and  K.B.M.Q.  Zaman,  "An  experimental  study  of  organized  motion  in  the  turbulent  plane  mix¬ 
ing  layer,"  J.  Fluid  Mech.  159  (1985),  85-104. 
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